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Fluorescent imaging plays a critical role in a myriad of scientific endeavors, particularly in the 
biological sciences. Three-dimensional imaging of fluorescent intensity often requires serial data 
acquisition, that is voxel-by-voxel collection of fluorescent light emitted throughout the specimen 
with a non-imaging single-element detector. While non-imaging fluorescence detection offers some 
measure of scattering robustness, the rate at which dynamic specimens can be imaged is severely 
limited. Other fluorescent imaging techniques utilize imaging detection to enhance collection rates. 

A notable example is light-sheet fluorescence microscopy, also known as selective-plane illumina¬ 
tion microscopy (SPIM), which illuminates a large region within the specimen and collects emitted 
fluorescent light at an angle either perpendicular or oblique to the illumination light sheet. Unfortu¬ 
nately, scattering of the emitted fluorescent light can cause blurring of the collected images in highly 
turbid biological media. We recently introduced an imaging technique called coherent holographic 
image reconstruction by phase transfer (CHIRPT) that combines light-sheet-like illumination with 
non-imaging fluorescent light detection. By combining the speed of light-sheet illumination with 
the scattering robustness of non-imaging detection, CHIRPT is poised to have a dramatic impact 
on biological imaging, particularly for in vivo preparations. Here we present the mathematical for¬ 
malism for CHIRPT imaging under spatially coherent illumination and present experimental data 
that verihes the theoretical model. 


I. BACKGROUND 


Image formation by temporal modulation of illumination light has emerged as a promising method to circumvent 
current limitations in biological imaging applications CHII]. In particular, image acquisition speeds can be greatly 
enhanced by overcoming the need to serially collect three-dimensional image data, as is required for techniques such as 
confocal and multiphoton laser scanning microscopy (MPLSM) |121113| . Imaging systems relying on temporal modu¬ 
lations of the illuminating light to form images collect image data from multiple points in the specimen simultaneously 
with a single-element detector. 

Recently, we introduced an imaging method we call CHIRPT, for coherent holographic image reconstruction by 
phase transfer, which allows coherent imaging methods, including holographic propagation, to be directly applied 
to incoherent contrast mechanisms such as fluorescence m- CHIRPT has demonstrated an ability to dramatically 
increase imaging speeds and the depth of held when compared to conventional biological imaging methods. 

In this Paper, we rigorously investigate the process of image formation in CHIRPT by an angular spectrum rep¬ 
resentation m to compute to the illumination intensity pattern. We then compute the transfer function and point 
spread function (PSF) of an idealized CHIRPT imaging system, and show how this allows the spatial phase of the 
imaging system to be recovered with a non-iterative algorithm. Finally, we verify the theoretical model with an 
experimentally determined transfer function. 
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II. PRINCIPLES OF CHIRPT IMAGING 

CHIRPT forms images by illuminating a large area in the object region with a temporally modulated intensity 
pattern that is unique to each point in the lateral-axial plane, (x, z). The contrast mechanism, which can be coherent 
or incoherent, responds to this modulation pattern and re-emits, transmits, or scatters signal light that has the same 
spatio-temporal structure as the illumination intensity. The spatio-temporally modulated signal light is measured 
with a single-element photodetector, such as a photomultiplier tube (PMT), providing a periodic voltage signal, S{t). 
The electronic spectrum of this signal, S{v), where v is the temporal modulation frequency, encodes a one-dimensional 
image of the object while preserving the spatial phase. A two-dimensional image in the (x, z) plane is reconstructed 
by coherently propagating the contrast intensity, which relies on the spatial phase of the illumination fields preserved 
in S{t). 

A CHIRPT microscope is represented schematically in Fig. A collimated, spatially-coherent illumination beam is 
brought to a line focus on the face of a circular modulation mask using a cylindrical lens of focal length Fc^y, oriented 
to focus the beam in the vertical dimension (y). The mask is designed such that the modulation frequency imparted 
to the line focus varies linearly with a function of radial position, R (Fig. &)• The varying spatial frequency presented 
to the line focus, varies as the mask is rotated at a constant angular velocity, causing the illumination beam 

to be diffracted from the mask with an angle that varies with time, Oi{t) = sin“^ [X fx^i{t)/ni]. An imaging system is 
used to image relay the modulation plane to the specimen. Schematically, we represent the imaging system as a tube 
lens of focal length Ft and an objective lens of focal length Fo, providing a system magnification of M = Ft/Fo- 



FIG. 1: Schematic representation of a CHIRPT microscope, (a) Image relay system in the (x, z) plane, (b) A 
modulation mask of the form m{R, tp) = 1/2 -|- (1/2) sgn [cos(AA: with Ak = 3/mm. (c) To provide a unique 

mapping of each (x, z) location in the object region, only the undiffracted beam and one of the first-order diffracted 
beams are permitted to propagate to the object plane. A horizontally-oriented slit is placed in the pupil plane to 

completely block one of the diffracted beams. 


The critical component of a CHIRPT imaging system is the illumination intensity pattern, /iii(r, t), which illuminates 
every position in the (x, z) plane in the specimen with a unique intensity modulation as a function the rotation angle 
of the mask, foit)- For the illumination intensity in the object region to meet this criteria, only the undiffracted beam 
and one of the first-order diffracted beams are permitted to interfere in the specimen region. This is accomplished 
with a horizontally-oriented slit placed in the conjugate image plane, which is possible due to a constant non-zero 
spatial frequency imparted to the illumination beam in the vertical dimension, y (Fig. &)• This spatial frequency is 
computed in Section [HI B| 

As the mask completes a full rotation, the full numerical aperture of the imaging system is scanned by the diffracted 
beam. Each temporal measurement on the photodetector corresponds to a specific angle of the diffracted beam in 
the object region, 02 (t)- The diffracted beam interferes with the undiffracted beam and creates an interference patten 
with a single spatial frequency in the lateral dimension, This produces the temporal photodetector signal as 

a function of time corresponding to a measurement of the lateral spatial frequency content of the specimen, and a 
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one-dimensional image in the lateral dimension can be recovered by a Fourier transform. As we show in the following 
analysis, spatial phase information is also encoded into this signal, and a two-dimensional image in the (x, z) plane 
can be recovered from a temporal measurement that corresponds to a single rotation of the modulation mask. 


III. THEORETICAL ANALYSIS 


CHIRPT encodes complex spatial image information into a temporal intensity by measuring the radiant flux from 
the contrast distribution on a single-element photodetector as the illumination intensity pattern is modulated over 
time. The CHIRPT signal has the form of a projection: 


S(t) = 



/iu(r,t) C(r,t) 




( 1 ) 


where t) is the illumination pattern in the object region, and C(r,t) is the contrast distribution. The contrast 
distribution describes the method by which illumination intensity is transferred to contrast intensity. For example, 
the contrast distribution can simply represent the transmittance of the object if the measured contrast intensity is 
the illumination source after propagation through an object displaying absorption. In that case, C{r,t) = T(r,t) = 
1 — A(r, t), where T is the intensity transmission of the object, and A is the intensity absorption of the object, and an 
image is formed from the object transmission. The contrast can also represent the spatial concentration distribution 
of fluorophores in an object, scaled by the absorption cross section for the excitation light and efficiency of fluorescent 
emission. Throughout this analysis, we shall assume that the contrast function is constant for the duration of a 
scan, so that C{r,t) —>■ C(r). We also assume that the contrast function is linearly proportional to the illumination 
intensity. 

The propagation phase imparted to the intensity emerging from the object arises from a change in the instantaneous 
modulation frequency at a given point in the contrast distribution as a function of scan time. In general, the spatial 
phase of the illumination source is encoded in the signal from the photodetector as a phase variation on the temporal 
carrier frequency. 

This phase transfer is a critical aspect of CHIRPT. Here, we develop a theory that describes the generation of 
CHIRPT signals. To derive the form of the signal generated by the photodetector, we consider a spatially-coherent 
line focus passing through a spatio-temporally varying modulation mask. Diffraction by the mask produces transverse 
wave vectors that vary with time. We model this diffraction by a set of plane waves propagating at angles with respect 
to the optic axis (z) that vary with time (Fig.jll. Using these angles, an expression for the electric held diffracted from 
the modulation mask for an arbitrary set of diffracted fields is derived. The phase accumulated upon propagation 
into the object region is computed via angular spectrum formalism applied to the re-imaged diffracted waves that 
interfere to produce the illumination intensity, taking the system magnification into account. An expression for the 
total illumination held in the object region is found by taking the square-magnitude of the sum of the diffracted 
orders. A pupil phase is included to account for systematic aberrations to the imaging system. 


A. Plane Wave Expansion of the Illumination Intensity Pattern 

CHIRPT is part of a broader class of imaging in which multiple diffractive orders are used to form a time-dependent 
illumination intensity pattern which encodes image information into a whole-held measurement. Other methods that 
fall in this category include synthetic aperture microscopy (SAM) jj] ITSUIS] and spatial-frequency-modulation for 
imaging (SPIFI) [5HZ]- As such, we begin by deriving general expressions for the illumination intensity in the object 
region when a spatially coherent illumination source is modulated by an arbitrary mask pattern. 

In the following analysis, we use the subscript pair (j, s) to denote the diffracted order and either the mask 
region with s = 1 or the object region with s = 2. For example, the set of k-vectors corresponding to the diffracted 
beams in the mask region are described by wavevectors kj i = ® 2 - Spatial coordinates are 

described with an arbitrary position vector, r = xe^ + yey + zBz- The optic axis of the system is set to be the 
z-direction (e^), and we define transverse k-vectors (in the x,y plane) by the notation 
and the position in the transverse plane as rj_ = x -I- y ej,. 

For simplicity, we assume a uniform plane wave of unit amplitude incident on a diffraction mask, and propagating 
along the optic axis, which is normal to the plane that defines the diffraction mask. The incident held contains a 
single spatial frequency, so the spatial portion of the electric held incident on the modulation mask can be written as: 
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i;inc(r)=e‘‘^^'- = e''=^^ (2) 

where ki = 2-^711 jX is the wavenumber of the incident wave and the wave vector is strictly along z: ki = kie^. 


The modulation mask causes the incident field to diffract into multiple directions propagating at various angles with 
respect to the optic axis. The structure of the mask dictates the direction of the diffracted wave vectors by imparting 
a transverse wavevector to the incident plane wave. Expressing the magnitude of the k-vector in the transverse plane 
{x,y) of a diffracted order j as = k±j^s(t), the magnitude of the axial k-vector for each diffracted wave in 

the mask region is: 


= |kzj.i(i)| = •klj-i(t) = kiJl- 


ki 


1 2 


We can therefore write the electric field for the j**' diffracted wave in the mask region as: 

Vj^i{r,t) = Qj = Qj exp [ikj_j- 4 (t) • rj_] exp [ikzj^i{t) z ]. 


( 3 ) 


( 4 ) 


where the amplitude of the diffracted field can be a complex quantity, i.e., aj = |aj |e'^“j . Physically, the phase of 


— lUj 

the amplitude, Za,, corresponds to a lateral translation of the modulation mask. 


The electric field distribution directly behind the modulator is image relayed to the object plane with an imaging 
system of magnification M (Fig.[^. Therefore the axial wavevector in the object plane is: 


kz,j,2{t) = \lkl- [fc_L,y2(t)]^ = k2 \ I- 


712 M k±^j^i{t) 
ni k2 


( 5 ) 


where we have used the fact that kj_,j^2{t) = ( 712 /^ 11 ) M /c_L.yi(t), and the wavenumber in the object region is scaled 
by the index of the immersion medium, 772 . Thus the electric field for the j**' diffracted wave in the object plane 
becomes: 


Vj ,2 (r, t) = Uj exp 


i — Mk_Lj,i(t) • rj^ 

, Til 


exp 



772 M k±j^i{t) 
771 /C2 


( 6 ) 


While Eq. ([^ accounts for the phase caused by diffraction away from the focal plane (system defocus), there is, in 
general, an arbitrary pupil phase, 4>pupii [k_Lj (t)], that must be included to account for other spatial phase variations 
imparted by the imaging system. Including the pupil phase, the expression for a diffracted order in the object region 
is: 


Vj^ 2 {'r,t) = aj exp [ik±j^ 2 {t) ■ r_L] exp [iA:^j, 2 (<) z] exp{i$pupii [k_Lj(t)]} 


( 7 ) 


In the object region, the diffracted beams {j = ±1) propagate at angles ±02(i) = sin“^ [±AA:j_j_2(0/(27’‘)], so that 
we need to account for diffracted waves at ±kx,j, 2 (i). A given mask can generate many diffracted waves, which can 
either be derived from higher diffracted orders of a particular mask lattice vector, or from the inclusion of multiple 
lattice vectors that can diffract beams to additional angles m- If we suppose that we have a total of N diffracted 
waves of appreciable power, then we have a total of 2N + 1 waves (including the zero-order term) that must be 
accounted for. The total illumination field propagating in the object region can thus be compactly expressed as 

N 

Eiu{r,t)= ^'i.2(r,^) (8) 

j=-N 


Since the mask is real and positive semidefinite, all of the diffracted waves for given index j have conjugate 
amplitude, i.e., —aj = a*. Also, each pair of diffracted waves ±j have the same axial wave vector, but conjugate 
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transverse wavevectors, so conjugate pairs in the summation can be paired: 
Siu(r,t) =aoe‘'==^e‘'*’p"p“(°) 


gi ^pupii[kx,j (t)] 

Here the j = 0 term is the undiffracted beam. The illumination intensity in the object region is the square modulus 
of the expression above, t) = t)|^. Defining the complex conjugate of the electric field as: 

Uj,2{r,t) = [vj^2{r,t)]* ( 10 ) 

the illumination intensity can be expressed as: 


N 

+2^ I Oj I cos 
1=1 


— M kj_ 1 1 (t) • r + Zaj 
m 


exp< ik 2 Z \ 1 — 


712 

Til k2 


1 2 


N N 

lill ^7.2 (r, t) 77fe,2 (r, t) 

j^-N k=^-N 


( 11 ) 


where N is the highest diffracted order included in the intensity pattern. 


B. Plane Wave Expansion of the Modulation Mask 


Utilizing Eq. ([^, the illumination intensity in the object region for an arbitrary modulation mask can be found 
by computing the set of transverse wavevectors kj_j_i(t) imparted on the illumination beam by the mask. To obtain 
an expression for the transverse wavevectors, we consider the mask as a decomposition of plane wave vectors. We 
assume that the mask is constructed of periodic modulations along directions described by lattice vectors Bp. We 
further assume that the mask modulates the amplitude of the field, requiring that it be positive semidefinite. With 
both of these assumptions in hand we can write the general form of the mask as a Fourier cosine decomposition of 
wavevectors: 


m(r, t) = bo + y] bp cos [kp(t) • r+ ap(t)] (12) 

23=1 

where we have utilized the fact that kp(t) = kp{t) Bp, and ap(t) is a phase shift that accounts for lateral translations 
of the modulation mask as a function of scan time. The term bo is simply a constant offset that must be selected to 
allow the mask to be positive semidefinite. All other terms correspond to periodic patterns in the modulation mask. 

Since the modulation mask is confined to the lateral plane, i.e., the mask provides no axial modulations, this 
wavevector is equivalent to the transverse wavevector experienced by the incident electric field. Therefore, kp(t) = 
kx,j,i(t), and we arrive at a general expression for the modulation mask in terms of the lateral wavevectors imparted 
onto the incident coherent illumination source: 

N 

m{r±,t) =bo + '^ bj cos[kxj,i(t) ■ r± + aj{t)] (13) 

i=i 


For j ^ 0, bj = 2\aj\. 


IV. CHIRPT IMAGING 


We now consider the form of the modulation mask used for CHIRPT imaging. In polar coordinates, {R,(p), the 
pattern of the CHIRPT modulation mask has the form: 

™(^3 = \ + isgn[cos(Afc R ^p)\ (14) 

where sgn[-] is included to account for the binary amplitude modulation resulting from the printing process, and A/c 
sets the highest density of the printed pattern. An example mask is shown in Fig. [Tj7. 
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We wish to decompose the transverse wave vector imparted on the cylindrically-focused beam into Cartesian 
components, where x is the lateral dimension, parallel to the line focus, and y is the vertical direction, perpendicular 
to the line focus. Equation Eq. (14) represents the full two-dimensional pattern of the mask, while we are interested 


in the region sampled by the line focus, located at j/ = 0, and approximately centered on one side of the mask such 
that the illumination distribution in x lies between the inner and outer radii of the mask (Fig. &)• To account for 
the changing modulation pattern sampled by the line focus during mask rotation, and thus find the local spatial 
frequencies as a function of scan time, we consider the mask: 


m[R,y) - ^po{t)\ = i -k ^sgn(cos{Afci? [(/? - v3o(i)]}) 


(15) 


where <po(i) describes the rotation angle of the mask. Since the mask rotates at a constant angular velocity, ipo{t) = 
2TTi'rt, where Vr is the rotational frequency of the mask. 

To uniquely determine the spatial phase of the illumination microscope, only the first diffracted order from the 
mask is permitted to interfere with the undiffracted beam in the object plane. Higher diffracted orders resulting from 
the binary modulation scheme are omitted in the following analysis, and also excluded in our CHIRPT experimental 
design, so that the mask we consider is: 

m[R,ip- ipo{t)] = ^ ^ cos{Afc R[tp- ipo{t)]} 


= ^ + ^cos[(j){R,ip,t)] 


(16) 


The local linear spatial frequency imparted by the mask in Cartesian coordinates is m-- 

^ dx 

t us 1 0 

^ dy 


(17) 

(18) 


The region sampled by the beam is along the line where = 0, so the local angular spatial frequencies can be 
calculated as: 


^X,l(^) - 

ky,l{t) = 


d(l}(R,(p,t) 1 . dd>{R,(p,t) 

cosip ---sm(/? 


dr r ^ dip 

. d(j){R,p,t) 1 d(j){R,p,t) 

SlUif ---h - COSif -p;- 

or r o ip 


- ^p—O 


(19) 

( 20 ) 


- p—O 


where we have converted linear spatial frequency to angular spatial frequency by multiplying by 27r. The local spatial 
frequencies evaluate to: 


kx, 1 (t) = 2 tt Akvrt = 2Tr Kit ( 21 ) 

ky, i(t) = Ak (22) 

and the transverse wave vector for the fundamental mask frequency is kj^_i^i(t) = 27r ki t e^, -|- Akey. The parameter 
Ki is the so-called chirp parameter, which relates the modulation frequency to lateral position on the modulation 
mask. 


The expressions above elucidate the key elements of CHIRPT imaging. In the lateral dimension, spatio-temporal 
modulations encode the lateral position with unique modulation frequencies. Moreover, each measurement in time 
corresponds to a unique lateral spatial frequency, demonstrating how CHIRPT measures the lateral spatial frequency 
content of the object. 


Returning to the general form of the modulation mask in Eq. (13), the transverse wavevector can be decomposed 
into Cartesian components to give: 


N 

m(x, y, t) = bn + bj cos[j 2tt Kitx + j Ak y Oj (t)] (23) 

i=i 

Since CHIRPT only utilizes one of the first-order diffracted modes and the undiffracted beam, the plane wave expansion 
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of the mask can be simplified further: 

m{x, y, t) = bo + bi cospTr mix + Ak y + ai (t)] (24) 

Since the mask contains a finite set of spatial frequencies, the function describing the mask must be constrained to a 
finite time window. Thus we can write the mask as: 


m{x,y,t) = w{t) {bo + bi cos[27rKi tx + Aky + ai(t)]} 


(25) 


where the temporal window function w{t) restricts the mask function to a single rotation of the disk and is used 
to describe apodization of the throughput energy. Throughout this analysis we assume that w{t) is a real-valued 
function. 

As noted previously, the CHIRPT illumination intensity pattern is formed by allowing only one of the diffracted 
orders to interfere with the undiffracted beam in the object region. This is accomplished with a horizontally-oriented 
slit placed near the pupil plane of the illumination lens (Fig. that passes only a portion of the j = 0 and 
j = 1 beams, which is possible due to a non-zero spatial frequency component in the vertical (e^) direction that is 
independent of the rotation angle of the mask [Eq. (22l]. Including only the j = -1-1 and j = 0 orders the illumination 
field [Eq. 


is F;iu(r, t) = no, 2 ( 1 ', t) + ni, 2 (r, t), where: 

^^ 0 , 2 ( 1 ',*) = oo exp [i/c 2 z] exp [i 4'pupii(0)] 

m 2 (r, t) = ai w{t) exp ( i27r — Mki t a; ) exp ( i — M Aky 
\ ni J \ ni 


(26) 


X exp 


^ \ (n2 MAk^ ^ 

lk2Z\ll-{ - MKit] -- - - 

m / Vni k2 


(27) 


X exp [i4>pupii(t)] exp [iai(t)] 

Note that w{t) does not apply to no, 2 (i’,t) because this field is stationary with scan time. 

As indicated in Fig. the optic axis of the illumination beam is not collinear with the center of the mask 
modulation pattern. The consequence is that the illumination intensity, and hence the CHIRPT signal, contains a 
non-zero carrier frequency. This is a critical aspect of CHIRPT imaging, since a non-zero carrier frequency permits 
separation of the positive and negative sidebands in modulation frequency space, allowing single-shot data acquisition. 
Following the notation outlined in Fig.[^, we let the lateral position in Eq. Eq. (271 be a: = Xc + Ax, where Xc denotes 
the location of the centroid of the illumination beam intensity in the lateral dimension in the mask plane. Defining 
the carrier frequency as r'c = in 2 /ni) M ki Xc, we can rewrite Eq. Eq. (27l as: 


vi 2 (r, t) = oi w{t) exp (i 27r t) exp ( i 27r — M ki t Aa: ) exp (i — M Ak y 

' m J \ ni 


, / A _ fn2 MAk 

X exp ik 2 z \ 1 — [ — M Kit ] —- - - 

V V^i / V”-i ^2 

X exp [i$pupii(<)] exp [iai(t)] 


(28) 


The illumination intensity is defined as the square modulus of the illumination field, /iii(r, t) = |£’in(r, t)\^. We may 
therefore write the illumination intensity in the object region as: 

/iii(r,t) = vo, 2 {r,t)uo, 2 {'r,t) + ni, 2 (r, t)ui, 2 (r, t) -k vo, 2 {r,t)ui^ 2 {r,t) + uo, 2 {r,t)via{v,t) (29) 

For brevity, we write the electric field for each beam as: 

no,2(r,t)=aoe‘^‘>.= (--’‘) (30) 

ni,2(r,t) =aiu;(f)e‘^i-^("’*) (31) 

where the overall phase argument of diffracted beam j is 4)^,2 (r,t). The illumination intensity becomes: 

/iii(r,t) = ttg -k a\w‘^{t) -k 2aoai w(t) cos [4>i,2(r,t) - 4>o,2(r,t)] 

= 0,0 + a\w‘^{t) -k 2 no oi w(t) cos [A$(r, t)] 


( 32 ) 













Equation Eq. (321 elucidates the key principle of CHIRPT imaging - namely, that the illumination intensity, and 
hence the CHIRPT signal, is formed by the sum to two separate intensity components: one with no dependence on 
the phase of the illumination beams, and therefore no carrier frequency, and another with spatio-temporal intensity 
fringes that depend on the phase difference between the two illumination beams, A$(r,t). Since the dependence 
of the CHIRPT signal on the phase difference of the illumination beams is independent of the contrast mechanism, 
coherent imaging methods can be applied directly to incoherent light for the first time. 

The overall phase of each beam results from the varied optical path length (OPL) that each beam experiences while 
propagating from the modulation mask to a position r in the object region. The OPL and the phase are related by 
the wavelength of the illumination light by: 


$,- 2 (r,t) = y OPL,- 2 (r,t) (33) 

Clearly the intensity modulation pattern is driven by the difference in OPL experienced by the two illumination 
beams: 


O-TT 27 r 

A$(r,t) = y [OPLi, 2 (r,t) - OPLo,2(r,0] = y^OPL(r,f) 


(34) 


Using the phase variations for each beam defined above, we can write the phase difference encoded into the illumination 
intensity as: 

A$(r, t) = 2Tri'ct + 2 tt — M Kit Ax H— - M Ak y 

rii Til 


+k 2 z 


' 1 - 


A 

ni 


M Kit] — 


712 M Ak 

ni /c2 


- 1 


+ A<I>o(t) -|- ai{t) 


(35) 


where A<i)a(t) = $pupii(27r ki t) — 4>pupii(0) is the aberration phase, the pupil phase at time t referenced to the on-axis 
pupil phase value. 

The illumination system is designed to allow portions of the j = 0 and j = 1 beams to pass at the equator of the 
objective lens (vertically), thereby removing the shear between the diffracted and undiffracted beams in the vertical 
dimension. As such, one can make the approximation that in the object region, Ak —^ 0 in the y-direction only, and 
the phase difference no longer varies with respect to the vertical dimension: 


A$(a;, z,t) = 2TTi'ct + 2 tt — M Kit Ax + k 2 z 

n-i 



1 -|- A<i>a(t) -|-Q!l(t) 


(36) 


Since the carrier frequency Kc only appears in the phase difference, we can split the illumination intensity into 
a DC term and an AC term, which we denote as Io{t) and Ii{r,t) respectively. The illumination intensity is then 
/in(r,<) = Io{t) +Ii{r,t), where: 

Io{t) = wo.2(r, t) uo,2{r, t) + ui,2(r, t) 7ti,2(r, t) 

= a^ +a\w‘^{t) (37) 

and 


/i(r, t) = vo, 2 {r, t) ui, 2 (r, t) -k vo, 2 {r, t) Mi, 2 (r, t) 

= 2 oo oi w{t) cos [A$(r, t)] (38) 

Note that the temporal dependence in Iq (t) is due only to the temporal window, which accounts for a loss in transmitted 
amplitude of the j = 1 beam with scan time. 

Although the mathematical formalism at this point corresponds to a sinusoidal amplitude grating, i.e., N = 1, 
experimentally the modulation mask is binary. We therefore select the amplitude coefficients to reflect the diffractive 
amplitudes from a square grating: oq = 1/2, and oi = I/tt. Thus the illumination intensities are: 

loit) = I + ^w'^{t) 

4 TT^ 

/i(r,t) = —w{t) cos [A$(r, t)] 

TT 


(39) 

(40) 
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We can also express the measured signal as a sum of the AC and DC components: 

S{t) = + hir,t)] C(r))^ = So{t) + Si{t) 


(41) 


The AC component of the illumination intensity is responsible for encoding the image into the temporal measurement 
from the single-element detector: 


Si{t) = (/i(r,t)C'(r))r = -w{t) (cos[A$(r,t)] C'(r))^ 

TT 


(42) 


The illumination intensity pattern driven by the phase difference in Eq. (351 varies with scan time in x and z, 


but is constant with respect to the vertical coordinate, y. This is expected, as the vertical spatial frequency on the 
modulation mask is constant with respect to time [Eq. ( |22| )] - a result of being invariant with mask rotation angle, 
<po(t)- The effect of the vertical spatial frequency is to impart a linear phase delay in the vertical dimension that 
is directly proportional to the mask density, Afc, and another linear phase delay that arises from the shear of the 
re-imaged diffracted order with respect to the optic axis. The linear phase delay with respect to the axial coordinate 
leads to an intensity modulation in the case where the full y spatial frequency extent of the zero-order and first 
diffracted order beam are allowed to pass through the CHIRPT illumination optical system. Placement of the filter 
that selects a portion of the y spatial frequency content eliminates this modulation that arises from the vertical spatial 
frequency shear (Fig. Ic). The presence of the restrictive spatial filter leads to the illumination intensity distribution 
given in Eq. (401 with a phase difference defined by Eq. (36 1 . 


Finally, since only the lateral and axial components of the contrast function are encoded into the temporal signal, 
we assume that the contrast function is separable with respect to the vertical dimension, i.e., C(r) = Cy{y) C{x, z), 
so that the AC signal is: 


-| pOO p pOO 

Si{t) = -w(t) / dy Cy{y) // dzdxC{x,z) 

J —oo J J —oo 


X cos { 2tt Vet + ^TT ( — M Kit ] Ax + k 2 z 


'I - { —MKit] - 1 


+ A4>tj(t) -I- oii{t) 


(43) 


Integration over the vertical coordinate in the expression above results in a constant multiplicative factor. Let us 
define an overall prefactor to account for this integration and the leading factor of I/tt as 7 = (I/tt) f dyCy(y). 

The signal collected from the photodetector must be real-valued, and thus it follows that the signal displays 
conjugate spectral symmetry. Both the positive and negative fundamental modulation sidebands, corresponding to 
the AC component of S(t), encode phase information, although they are conjugates to one another. For image 
reconstruction, only one such signal is needed. The complex CHIRPT spatial frequency spectrum from the temporal 
data trace is isolated by filtering out the positive frequency sideband of the measured signal, then inverse Fourier 
transformed. 

After these operations, we arrive at an expression for the measured spatial frequency distribution; the amplitude of 
the the time-dependent trace contains the object intensity spatial frequency distribution, and the phase represents the 
propagation and modulation phase of the illumination intensity accumulated between the modulator and the object 
region. In the analytic theory presented here, this is readily observed by representing the cosine term in Eq. (431 as 


a sum of complex exponentials. To simplify the notation in the following expressions, we define the chirp parameter 
in the object region as: K2 = (n2/ni) M ki, and observe that the lateral spatial frequency in the object region is 
fx, 2 it) = K 2 t = kx, 2 {t )/Then we can write: 


Si{t) = -w{t)^ jj dz dx C{x, z) i2TT Vet + ikx, 2 (t) Ax 


+ ik2Z 


A- 


kx,2{i) 


- 1 + i A$a(f) +iai{t) 


-I- c.c. 


(44) 


where c.c. denotes the complex conjugate. We represent the positive and negative temporal sidebands centered at 
±Vc as Si+{t) and Si-{t) respectively, such that Si{t) = 5'i+(t) -I- Si-{t). Each term results from a single sideband 
in the modulation frequency domain, and is therefore a complex quantity. 

While the intensity is by definition real-valued, each sideband can be thought of as resulting from the individual 
complex components of the intensity. Denoting the complex components of the illumination intensity in the same 
manner as the measured signal, i.e., /i(r,t) = /i+(r,t) -|-/i_(r,t), where Ii{r,t) G K and /i±(r,t) G C, we can write: 
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Si±{t) = (/i±(r,t)C(r))^ = {vi, 2 {r,t) uo, 2 (r, t) C'(r))^ 


(45) 


Notice that in order for this statement to hold, the complex components of the intensity must be conjugate to one 
another, i.e., /i 4 .(r,t) = [/i_(r,t)]*, and consequently the positive and negative temporal sidebands carry redundant 
spatial phase information. In practice, we select the positive temporal first harmonic for image reconstruction. Finally, 
we arrive at an expression for the image from a single lateral-axial plane encoded by CHIRPT: 




JJ dz dx C{x, z) exp}.ikx, 2 {t) + ik 2 z 


' 1 - 


kx,2{t) 


- 1 


(46) 


A. Spatial frequency projections 


CHIRPT encodes lateral image information by sequentially projecting an illumination pattern with a single lateral 
spatial frequency, kx, 2 (t), onto the specimen. The contrast function, C{x, z) is measured by collecting the total radiant 
flux emanating from the specimen region for each spatial frequency. The CHIRPT signal in time, Si+{t), represents 
a projection of the spatial frequency of the illumination intensity onto the object. Thus the CHIRPT signal can be 
recast into the lateral spatial frequency domain by relating it to the scan time. This is clear from Eq. . Rewriting 
the expression for the CHIRPT signal in lateral spatial frequency coordinates, we And: 


Si+{fx) = 


dz exp i fc 2 z 


' 1 - 


fx,2{t) 

f2 


1 2 


- 1 


/ oo 

da; C{x,z)e^^^f-^ 

-OO 


(47) 


where we have applied the shorthand notation fx = fx, 2 (t) to represent the lateral spatial frequency of the collected 
signal. Notice that in this representation, the temporal window function limits the spatial frequency content of the 
measured signal, which ultimately sets a limit on the numerical aperture of the measurement. 

The integration over the lateral coordinate has the form of a spatial Fourier transform of the contrast function with 
respect to the lateral coordinate, so the preceding expression can be written: 


Si+{fx) = i «;(/.) e‘giA$„(A)giai(A) 


dz exp i fc 2 z 


' 1 - 


fx,2{t) 

/2 


-1 H Cx{fx,z) 


(48) 


where Cx{fx^ x) is the one-dimensional Fourier transform of C{x, z) with respect to x. 

Equations Eq. and Eq. (481 are descriptive of CHIRPT imaging, in that they demonstrate how amplitude 


and phase information from the spatial frequency domain are encoded into a temporal measurement. By virtue of 
the phase delay with respect to the carrier frequency in the CHIRPT signal, the total spatial phase accumulated 
between the modulation plane and the object plane is encoded into a single temporal measurement. This includes 
phase variation from defocus, which is explicitly included in the theory presented above, as well as higher-order spatial 
frequency phase variations. 

Physically, one can think of the complex CHIRPT signal Si+it) ^ Si+{fx) as encoding the amplitude and phase of 
the lateral spatial frequency content, fx- A portion of the longitudinal spatial frequency content of the specimen is also 
encoded through the phase imparted by physical defocus. In other words, each temporal data point corresponds to the 
amplitude of a lateral spatial frequency fx, 2 {t), and a unique axial spatial frequency /z, 2(^)1 because the axial spatial 
frequency of a plane wave is determined by the transverse spatial frequency and the total wavenumber constrained 
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by the dispersion relationship dictated by the Helmholtz eqnation. This point is shown explicitly in 


Section IVT] 


V. POINT EMITTER SIGNAL 


It is instrnctive to stndy the case of a point emitter, which is formally eqnivalent to examining the transfer fnnction 
in CHIRPT imaging. The contrast fnnction for a point emitter located at (xp, Zp) in the (a;, z)-plane, relative to the 
centroid position (xc, Zq), is described by a prodnct of Dirac-(5 functions: C{x, z) = S{x — Xp)d{z — Zp). Note that for 
simplicity we have omitted proportionality constants such as the absorption cross section. Assuming no pupil phase 
for the imaging system (A<i)a(t) = 0) and zero translation of the mask with scan time (ai(t)), Eq. (46) reduces to: 


Si+{t) = - w{t) exp [i27r {i/c + K 2 Xp) t] exp 


'lk2 Zp 


' 1 - 


kx,2(t) 


- 1 


(49) 


It is clear that the temporal first harmonic of the photodetector signal encodes lateral position of a single emitter 
through carrier frequency, and axial position through frequency chirp. In the paraxial limit, where [kx^ 2 {t)/k 2 ] <C 1, 
we can make the approximation: 



1 - 


kx,2{t) 

k2 


= 1 - 



and Eq. (49) becomes: 


Si+{t) « - w(t) exp [i27r (vc + K 2 Xp) t] exp 


- ^ / \2 
—ITT - (K2 t) Zp 

n2 


(50) 


(51) 


Note that in the paraxial approximation, axial phase manifests in the form of a Fresnel zone plate (FZP) in the lateral 
spatial frequency domain, suggesting a strong analogy to off axis spatial holography Enna. 

For a collection of Np point emitters with locations given by a set of coordinates {{xp, Zp)}, the contrast function 
becomes: 


Np 

C{x, z) = GpS {x — Xp) S {z — Zp) (52) 

p=i 

and the signal from the photodetector is a linear combination of the signals from each emitter. In the paraxial 
approximation, the temporal scan encodes a summation of FZPs in the lateral spatial frequency domain because 
the total signal is a linear combination of the signal from each emitter |23| . Uniqueness of the illumination pattern 
for conjugate locations in the object plane, e.g., {xi,zi) and (xi,—zi), is maintained by the shift in instantaneous 
frequency relative to the carrier frequency, preventing ambiguity of emitter location in the measured signal. 


A. Instantaneous modulation frequency 


The temporal modulation frequency of the positive first harmonic signal is dependent upon the position of the point 
emitter in the lateral-axial plane. This can be made explicit by computing the form of the instantaneous modulation 
frequency, v{t)\ 


J A4>(xp,Zp,t) 
dt\ 2 tt 


d f AOPL{xp,Zp,t)\ 

1 A j 


= l/c + K2 Xp 


n2 



1 9A$„(t) 

27r dt 


1 dai{t) 
27r dt 


(53) 


This relationship between the position of the point emitter and the temporal modulation frequency of the measured 
signal is convenient for understanding how each position has a unique temporal modulation pattern in the object 
plane. 
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Figurej^shows the real portion of the complex temporal CHIRPT traces, Re computed for a single emitter 

with two different defocus values. For simplicity, we set the lateral position of the emitter to zero {xp = 0), assumed 
a perfect image relay system for which the aberration phase is zero, and set the phase due to lateral translations of 
the disk to zero. The signal was computed for a CHIRPT microscope with A = 532 nm, M = 95, Ak = 70/mm, 
and an air immersion illumination objective with NA = 0.8. While the typical carrier frequency for such a system 
would be on the order of 10-100 kHz (depending on the rotation frequency of the modulator mask), we set the carrier 
frequency to 2 kHz to make it possible to visualize the chirp in modulation frequency as a function of scan time. A 
Gaussian temporal window function with full width at 1/e point of 20 ms was also included. 


a 



Scan time (ms) 

FIG. 2: The instantaneous modulation frequency varies with the position of each emitter in the object region. 
Here we plot the real portion of the complex CHIRPT signal computed for (a) Zp = 4.9 /rm and (b) Zp = —1 fiui. 
Parameters of the imaging system simulated here can be found in the main text. 


The real component of the complex temporal signals in Fig. |^show how a one-dimensional temporal measurement 
encodes the emitter location in two dimensions by deviations from the carrier frequency. In the analytic example 
shown here, the magnitude and sign of frequency chirp determine the emitter location (xp, Zp). In general, other phase 
disturbances, such as systematic aberration phase, are also encoded into the modulation frequency chirp. In practice 
a single emitter can be used to empirically determine the systematic aberration phase of the microscope, allowing for 
digital aberration correction to be applied to images collected with CHRIPT without the need for iterative algorithms 

El- 

Each temporal CHIRPT trace in time represents a one-dimensional line image in the lateral dimension. The image 
is recovered by demodulating the temporal traces by the average frequency, converting the scan time to lateral spatial 
frequency via fx, 2 {t) = {n 2 /ni) M Kit, and computing the Fourier transform with respect to /x, 2 - Line images 
corresponding to temporal traces displayed in Fig. [^are shown in Fig. where a fast Fourier transform (FFT) was 
used to compute the images from the numeric temporal data. 


VI. CHIRPT TRANSFER FUNCTION 


The expression for the one-dimensional CHIRPT signal for a point emitter can be used to compute the two- 
dimensional transfer function in the spatial frequency domain, ifx,fz)- As shown in Eq. (471, the lateral spatial 
frequency, fx, can be substituted in place of the scan time, t. Making this substitution, we express Eq. (491 as: 


= 2 


i 27r fx Xp 


h Zp 



Sl+{fx) 


(54) 
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a 




FIG. 3: One-dimensional images computed from the complex CHIRPT data in Fig. Both the intensity of the 
contrast signal and the phase difference between the two illumination beams are encoded for each image, (a) The 
image computed with Zp = 4.9 /im shows positive concavity in the spatial phase, consistent with the positive chirp 
observed in Fig. (b) Similarly, the image computed for Zp = —1 fim shows negative concavity in the spatial 
phase. The width of the distribution is much more narrow in this image because the magnitude of the defocus is 

significantly smaller. 


where we have demodulated the signal by the carrier frequency for simplicity. The transfer function for CHIRPT 
imaging under single-frequency illumination can be computed with a spatial Fourier transform of this expression with 
respect to z. Since the lateral position of the emitter contributes only a linear phase ramp, we set Xp = 0 for simplicity; 
then the transfer function is: 


1 ro° ik2Z 

J dz e 

/*oo 


-1 


-J e 


—i 271 fz z 


1 r°° 

= ^ / dz exp[i27r (/^ - A/^) z] 

tX S {fz - ^fz) 


(55) 


where A/j, represents the difference in longitudinal spatial frequency between the diffracted and undiffracted beams, 
and (5(-) represents the Dirac-delta function. The Dirac-(5 function is non-zero only when its argument is equal to 
zero, so the transfer function is zero everywhere in the spatial frequency domain except for where fz = Afz ■ Thus we 
can write the longitudinal spatial frequency in the object domain, which is constrained by the dispersion relationship 
imposed by the Helmholtz equation for plane wave propagation: 


fzAt) = y 



[A , J 

2 1 

1 

1—1 

—fxAt) 
L«2 j 

1 —I 

I 


(56) 


Since there is a unique mapping of the lateral and axial spatial frequency information, each temporal measurement 
Si^{t) defines a two-dimensional plane wave. Summation of all the plane waves over the duration of a temporal scan 
forms a 2D image from a ID data trace. 


Equation Eq. (551 elucidates a fundamental limitation of CHIRPT imaging, and coherent diffractive imaging in 
general - the transfer function contains vanishingly small longitudinal frequency support, and thereby the imaging 
system does not display axial sectioning. In his seminal paper from 1969, E. Wolf notes that holography captures a 
portion of the complex object response that is determined by the illumination conditions |24| . Since that information 
is confined to the surface of the Ewald sphere in the spatial frequency domain, the possibility of three-dimensional 
image reconstruction with a single view is prohibited. The radius of the Ewald sphere is defined by the wavenumber 
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of the light. Noting that the wavenumber can be related to the linear spatial frequency by ^2 = 211/2 = 27r(n2/A), 
Eq. (56) can be rearranged to show: 


/I — fx,2 + ifz,2 + /2)^ 


(57) 


which represents a slice through the Ewald sphere at fy = 0 and shifted in longitudinal frequency by / 2 . The Ewald 
circle described by this expression is shifted by /2 in the axial dimension because CHIRPT encodes the difference 
between the diffracted and undiffracted beams. Equation Eq. ^ represents the full Ewald sphere for an electric field 
diffracted from the mask. 


VII. EXPERIMENTAL CHIRPT DATA 


To validate the theoretical analysis presented here, we imaged a 100 nm-diameter fluorescent nanodiamond (END; 
ND-400NV-100nm-10mL, Adamas Nanotechnologies) as a function of defocus. Since the END was sub-diffraction 
limited for the CHIRPT microscope, it acted as a point emitter and allowed for the transfer function of the micro¬ 
scope to be determined. A laser with 532 nm nominal wavelength (Sprout 5G, Lighthouse Photonics) was brought to 
a horizontal line focus with a 150 mm focal length cylindrical lens (ACY254-150-A-ML, ThorLabs) onto the spinning 
modulator. A modulation pattern with Ak = 70/mm was printed in aluminum onto a glass disk (Projection Tech¬ 
nologies) . The modulation plane was image relayed to the specimen region in two stages to allow unimpeded access to 
a conjugate of the objective lens pupil plane for filtering with the horizontal slit. The image relay system consisted of 
125 mm and a 100 mm focal length lenses (ACH-254-125-ML and ACH-254-100-ML, ThorLabs) for the first stage, and 
a 250 mm focal length lens (ACH-254-250-ML, ThorLabs) and 0.8 NA objective lens (Zeiss, N-Achroplan 50x/0.8 NA 
Pol) for the second stage. The overall system magnification was 95x. Fluorescence was collected in the epi direction 
by reflecting the illumination beams into the specimen with a dichroic beamsplitter (FF562-Di03-25x36, Semrock) 
and collecting fluorescence transmitted through the dichroic. Fluorescence emission was Altered with a bandpass Alter 
(ET645/75m, Chroma). Fluorescence intensity was measured with a photomultiplier tube (PMT; H7422P-40, Hama¬ 
matsu). The signal was amplified with a low-noise current amplifier (DHPCA-100, Femto) and digitized with a data 
acquisition card (PCI-6110, National Instruments). Specimen positioning was controlled with a three-dimensional 
motorized stage ((x,y): MS-2000; z: LS-50, Applied Scientific Instruments). Data collection and stage control were 
implemented in a custom C=/= application written in house. 

Figure]^ shows the average OTF computed from 20 images of the END as a function of defocus, Az. The overlaid 
red line represents the Ewald phase, A/)Ewaid(/x,/z), computed with the experimental parameters described above. 
The excellent agreement between the measured and computed transfer function verifies the theoretical analysis. 

Amplitude (dB) 
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FIG. 4: Experimental transfer function from a CHIRPT microscope, measured with a sub-diffraction limited 

END. The theoretical Ewald phase is plotted in red. 


VIII. EMITTERS OF FINITE EXTENT 

So far we have computed the CHIRPT signal for point emitters. However this is only valid in the case where a 
physical emitter is smaller than the diffraction-limited resolution of a CHIRPT microscope. To explore CHIRPT 
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imaging in the case of finite-extent objects, we consider a contrast function defined by a two-dimensional Gaussian 
distribution in the (x, z) plane: 


C{x, z) = P exp < — 


{x-Xpf {z-ZpY 


wt 


wt 


(58) 


Note that factors such as the absorption cross section are absorbed into the scalar constant, /?. Inserting this contrast 
function in Eq. (46), we find: 


Si+{t) = i w{t) ttPw^w^ e-“" e“" e^ (G. 2 (i)/fc 2 )" i] 


X e 


i 2tt Vc t gi (t) i 27r (k, 2 t) Xp ^ ^2 ij 


(59) 


The complex component of the temporal CHIRPT signal in Eq. (59) is identical to that of the signal for a point 
emitter we have previously derived [Eq. (49)]. In addition to this complex term, the amplitude of the temporal signal 
is modulated with a real function that depends upon the size of the emitter in both the lateral and axial dimensions, 
expressed through the Gaussian widths Wx and Wz respectively. As such we can write: 


^(Gaussian) = TT jS Wx Wz A{t; Wx, Wz) 


(60) 


where ^^^(t) is equivalent to the expression in Eq. (49), and the amplitude of the complex temporal signal depends 
on the parameters of the Gaussian emitter: 


A{t; Wx,Wz) = exp — (tt K 2 f)^ {wl — wf) 


exp { ^ (^2 Wz)'^ 


A- 


kxaX) 


-[ 2 


- 1 


(61) 


To understand how this amplitude modulation affects the measured image, it is useful to covert the notation to the 
lateral spatial frequency domain: 


A{fx-,Wx,Wz) = exp 





exp 



(62) 


In this form, the amplitude modulation describes the finite support of the measured GHIRPT image in the lateral 
spatial frequency domain. Several properties of this function are worth noting. 

First, there is a Gaussian dependence on the inverse of the width of the emitter in the lateral dimension. This is 
fully expected as the GHIRPT trace in time physically corresponds to a measurement of the lateral spatial frequency 
content of the object. As the width of the object increases, fewer lateral spatial frequencies are needed to reproduce 
an image of the object, and consequently the finite support is reduced. 

There is a similar dependence on the axial width of the object, although there is additional term that modifies the 
finite support which follows an Ewald-like form. The axial component of the finite support has a profound affect on 
images collected with GHIRPT, as it implies that the axial size of an object limits the lateral resolution of the image. 
For example, consider an object for which Wx —t 0, but Wz is larger than the diffraction limit. The frequency support 
becomes: 


lim A{fx-,Wx,Wz) 

Wp —>-0 


i (fes 


= e 



(63) 


Glearly the finite support is reduced by the amplitude modulations due to the finite width of the object in the axial 
dimension. 

Finally, we note that in the limit that both Wx —t 0 and Wz —t 0, the amplitude function, tends towards unity, 
and the temporal trace, when normalized by the object widths, reproduces the form computed for a point emitter 
previously. 


lim 


^(G, 

’1+ 


aussian 1 


it) = 


(64) 
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IX. COMPARISON WITH SPIFI 


We have shown that in order to achieve a unique mapping of axial location with lateral spatial phase only one 
diffracted beam is passed through the imaging system and allowed to interfere with the undiffracted beam in the 
object region. If both of the diffracted beams (j = ±1) are allowed to pass the imaging system, the CHIRPT 
microscope becomes equivalent to a SPIFI microscope EH?]. Due to the conjugate symmetry of the axial phase 
component for the diffracted orders, the transfer function for SPIFI is symmetric about the lateral spatial frequency 
axis. Consequently, the illumination pattern is not unique for all points in the {x, z) plane, and the axial position of 
an emitter cannot be uniquely determined. 

The temporal first harmonic of the photodetector signal for a point emitter measured with SPIFI can be expressed 
by considering the field combinations that contribute to the sideband centered at +i'c- The corresponding complex 
portion of the illumination intensity is /i+(r, t) = '^ 1 , 2 ( 1 *, i) ?^o, 2 (r, i) + ^«-i, 2 (r,f)^^o. 2 (r, t), and the complex SPIFI 
signal at the positive sideband is: 

Si+it) = (/i+(r, t) C{r))^ = (ui, 2 (r, t) uo, 2 (r, t) C'(r))^ + (M-i, 2 (r, t) uo, 2 (r, t) C{r))^ 

= ao ai w{t) (exp [i («'n 2 (i’, t) - 5>o,2(i’, t))] C'(r))^ 

+ao oi w{t) (exp [i (-$_i, 2 (r, t) + $o, 2 (r, t))] C(r))^ (65) 

Following the formalism derived in Section [rV| it is straightforward to show that the phase differences in the expression 
above have conjugate symmetry in both lateral spatial frequency and axial spatial frequency. In other words, instead 
of a single measurement in time corresponding to a unique lateral and axial frequency pair, {fx,fz)i the spatial 
frequency pair (/a,, —fz) is simultaneously encoded. Explicitly, the SPIFI signal for a point emitter is: 


Si+(t) = - w(t) exp [i27r (i/c + ^2 Xp) t] exp 


ik2Zp < W 1 - 


kx,2(t) 


1 2 


- 1 


+ - w{t) exp [i 27r (i/c + K 2 Xp) t] exp 


r 

' 1 -r 

-ik2Zp < 



kxAt) 


1 2 


- 1 


= w{t) exp [i 27r (t'c + K 2 Xp) t] cos 


k2Zp < \ 1- 


kx,2{t) 


n 2 


- 1 


( 66 ) 


From this expression it is clear that the uniqueness of the axial phase is destroyed if both diffracted orders are allowed 
to contribute to the illumination intensity pattern. As a result, a point emitter imaged with SPIFI results in a 
temporal pattern that is similar to the well known twin-image problem in holography, wherein an object is known to 
be defocused but the sign of the defocus is unknown. 

We computed the transfer function and PSF for CHIRPT and SPIFI imaging (Fig.|^. While the CHIRPT transfer 
function displays a unique mapping of each point in the spatial frequency domain, the SPIFI data shows that for each 
lateral spatial frequency there are two possible axial spatial frequencies of equal magnitude but opposite signs. This 
is a result of the cosine in the axial dimension. Similarly, the PSFs for CHIRPT and SPIFI are distinctly different - 
whereas CHIRPT has a smooth PSF about the axial dimension, the SPIFI PSF displays spatial modulations. 


X. DISCUSSION 


We have presented a theoretical analysis of CHRIPT imaging, a new form of optical microscopy that captures 
the spatial phase of a specimen with a measurement of incoherent (e.g., fluorescent) light. The analysis presented 
here assumes plane wave illumination with a single-frequency, spatially coherent source. The analytic expressions 
developed here are valid for CHIRPT microscope systems in which focusing in the vertical dimension can be ignored. 
We have verified the analytic expressions with experimental data. 

CHIRPT makes it possible to directly apply the powerful suite of coherent diffraction imaging (CDI) tools - which 
rely on spatial phase information available with coherent light - to images formed with incoherent light for the first 
time. CHIRPT opens a window to deep, large volume, and high speed in vivo imaging that will bring new scientific 
insights to a range of fields. These capabilities are poised to dramatically impact biological imaging in particular. 
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FIG. 5: Theoretical OTF and PSF for CHIRPT and SPIFF (a) The CHIRPT OTF follows an arc in the spatial 
frequency domain defined by the differential axial phase, Acjj-Ewaidifx, fz)- Here the numerical aperture in each 
dimension is defined by (b) Conversely, the SPIFI OTF contains two arcs with conjugate 

symmetry in longitudinal spatial frequency, (c) The CHIRPT PSF contains unique spatial phase information in the 
(x, z) plane, (d) The SPIFI PSF contains modulations in the intensity pattern which result from the degeneracy in 

axial spatial frequency content. 
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